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ABSTRACT 

We present a laboratory study of silicate glasses of astrophysically relevant 
compositions including olivines, pyroxenes and melilites. With emphasis on the 
classic Si-0 stretching feature near 10/zm, we compare infrared spectra of our new 
samples with laboratory spectra on ostensibly similar compositions, and also with 
synthetic silicate spectral data commonly used in dust modeling. Several different 
factors affect spectral features including sample chemistry (e.g., polymerization, 
Mg/Fe ratio, oxidation state and Al-content) and different sample preparation 
techniques lead to variations in porosity, density and water content. The convo- 
lution of chemical and physical effects makes it difficult to attribute changes in 
spectral parameters to any given variable. It is important that detailed chem- 
ical and structural characterization be provided along with laboratory spectra. 
In addition to composition and density, we measured the glass transition tem- 
peratures for the samples which place upper limits on the formation/processing 
temperatures of these solids in space. 

Popular synthetically-generated optical functions do not have spectral fea- 
tures that match any of our glass samples. However, the ~10/im feature gener- 
ated by the synthetic data rarely exactly matches the shape and peak position 
of astronomically observed silicate features. Our comparison with the synthetic 
spectra allows astronomers to determine likely candidates amongst our glass sam- 
ples for matching astronomical observations. 

Subject headings: stars: AGB and post-AGB — (stars:) circumstellar matter - 
(ISM:) dust, extinction — infrared: stars 
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Introduction 



Understanding the nature and formation of cosmic dust is crucial to our understand- 
ing of the cosmos. Over its 50-year history, infrared (IR) astronomy has shown that dust 
contributes to the physical processes inherent in star formation and mass-loss from evolved 
stars, as well as to several interstellar processes such as gas heating and the formatio n of 
molecules (e.g. IVideen fc Kocifajll2002l ; lDrainell2003l ; lKrishna-Swamyll2005t lKrugelll2008l ). In 
particular, silicate grains dominate dust emission in many astrophysical environments. The 
"amorphous" ~ 10/xm and ~ 18/im silicate spectral features have been observed in almost 
every direction and to almost any distance, but the precise nature of this silicate dust re- 
mains a mystery. Here we present a laboratory investigation of amorphous silicates, the type 
of dust grains most frequently inferred to exist from observational data. 



1.1. A brief history of the "amorphous" silicate spectral features 



The classic "10 /mi" silicate feature was first obse rved in the late sixt ies in the IR spectra 
of several M-type giants and red supergiants (RSGs; iGillett et al.lll968l ). Sho rtly thereafter 



a 10 / im absorption featur e was discovered in the interstellar medium (ISM; iKnacke et al. 



19691 ; lHackwell et al.lll970l ). Since then, it has been found to be almost ubiquitous, occurring 
in many as trophysical enviro nments including the solar system and extrasolar planetary sys- 



tems (e.g.. lMann et al.ll2006l . and references therein), the circumstellar regions of both young 



stellar objects and evolv ed intermediate mass stars (asymptotic giant branch; AGB stars, and 
planetary nebulae; e.g. JSpeck et al.ll2000l ; ICasassus et al.ll20011); m any lines of sight through 



the interstellar me dium in our own galaxy (e.g., IChiar et al. 



20071 ); and in nearby and dis- 



tant galaxies (e.g..lHao et al.ll2005t ). Initially this feature was attributed to silicate minerals 



( Woolf fc Ney 1969h , based on mixtures of spectra of crys talline silicate species predicted to 



form by theoretical models (jGaustadlll963 



Gilmanlll969l ). However, laboratory spectra of 



crystalline silicate minerals s howed more structure within the feat ure than observed in the 
astronomical spectra (see e.g.. lWooldll973l ; lHuffman &; St appll 19731). Subsequent c omparison 



with natural glasses (obsidian and basa l tic glass; from e.g., 



artificially disordered silicates ( IDaylll979l ; iKratschmer fc Huffman 



Po 



ack et al 

1979J) showed that 



19731 ) and with 
amor- 
phous" silicate was a better candidate for the 10 /im feature than any individual crystalline 
silicate mineral. Since then, it has been commonly assumed that description of silicate as 
disordered or "amorphous" is synonymous with glassy silicate. However, this is an oversim- 
plification. The term "glass" has a specific definition, i.e. the solid has no long-range order 
beyond nearest-neighbor atoms. "Crystalline" is often taken to mean single crystals, but it is 
possible to form poly-nano crystalline agglomerates, with a continuum which essentially ex- 
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tends from a true glass to a single crystal grain. Furthermore, natural and synthetic "glasses" 
often contain microliteg^In addition, one might expect agglomerated p articles to be po f ymin- 
eralic, and possibly to contain both crystalline and glassy constituents. iNuth &: Hechtl (119901 ) 
introduced the idea of "chaotic silicate" in which the level of disorder is even greater than 
for glass. A chaotic silicate does not have to be stoichiometric, can contain different compo- 
sitional zones within a single grain, and may be porous, and therefore much lower density 
than a glass. This range of possible grain types is demonstrated schematically in Fig. [TJ 

In astrophysical environments, whether a solid is glassy, crystalline or some combination 
of the two has implications for its formation, and subsequent processing, evolution and 
destruction, and thus it is important to have tools to distinguish between such grain types. 
For example, true glasses with no inclusions will not transform into crystals below their 
glass transition temperature (T ff ; see § ll.3p . whereas a glass already containing microlites 
can continue to crystallize at slightly lower temperatures. In the case of terrestrial obsidians, 
with 900 K < T g < 1000 K, ele mental diffusion pro files suggest that crystal growth can 



continue down to ~700K (e.g., IWatkins et al.l 120091 ) . If poly-nano-crystalline grains can 



be distinguished spectroscopically in the lab from t ruly glassy grains, we can test for their 
presence in astrophysical environments. For instance ISpeck fc Hofmeisterl (120041 ) showed that 
there is a difference between single crystal and polynano crystalline silicon carbide (SiC), 
while ISpeck. Thompson fc Hofmeisterl (120051 ) showed that glassy SiC looks different than 
various crystalline samples. These laboratory data were invoked to explain changes in SiC 
grains formed as carbon stars evolve. For silicates, laboratory data on poly-nano-crystalline 
samples are lacking. Whether grains are glassy or poly-nano-crystalline is an indicator of 
whether grains form or are processed above T g . If the two forms cannot be distinguished 
in laboratory spectra, then the "amorphous" nature of silicates in space would no longer 
necessarily imply a truly glassy structure, allowing the possibility of higher dust formation 
temperatures. 

Many observations have shown that the 10 fim silicat e feature varies from object to ob- 
ject and even within a single object both temporally (e.g.. lMonnier et al.lll998l ) and spatially 
(e.g. in rj Car, N. Smith, Pers. Comm). Within a single type of astrophysical object, the fea- 
ture sho ws huge variations in terms of its peak position, wid th and its ratio to the ~ 18/xm fea- 
ture (e.g. lSpecklll998l ; IOssenkopf. Henning. fc Mathislll992[ ). Variations in feature shape from 
star to star cannot be explained in terms of optical depth or grain temperature effects. Several 
interpretations of these observations have been suggested including: grain size effect (e.g., 



1 micro- or n ano- crystalline incl usions within a glassy matrix (see e.g., iPollack et al.l Il973t iJager et al 



1994J). Prior to IZachariasenl <| 1932T ) the difference between nano crystalline (i.e. ceramic) and glassy solids 



was not understood. 
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Papoular fe Pegouriall983l ); Mg/Fe ratio and ( Mg+Fe)/Si ratio (e.g. iDorschner et al.lll9 95); 



inclu s ion of oxide grains (e.g . JSpeck et aU2000f ); increasing crysta llinity (e.g. JSylvester et al. 



1998 



Bouwman et al.l 120011); grain shape (e.g.. iMin et al.l I2007D ; and grain porosity (e.g. 



Voshchinnikov fe Henning[2008 ; Henning fe Stognienko 1993 ). However, all models of these 



effects utilize laboratory spectra, and are only as reliable as the data that goes into them. 

The near ubiquity of "amorphous" silicate features and their variations in strength, 
width, peak position and the ratio of the 10/im/18/im features potentially provide the 
diagnostic tools to understand the detailed mechanisms by which dust is formed, processed 
and destroyed. However, existing laboratory and synthetic spectra are not sufficiently well 
understood to achieve this goal. 



1.2. A brief synopsis of existing laboratory and synthetic spectra for 

disordered silicates 



Since the discovery of the 10 /im feature, there have been many laboratory studies pro- 
ducing IR spectra and optical functions^ of various samples for comparison with and mod- 
eling of observational data. In addition, synthetic optical functions have been derived from 
observational spectr a, often combined with laboratory mineral data, in order to match the ob 



serve d features (e.g. , iDraine fe Leeill984l ; IVolk fe Kwoklll988l ; lOssenkopf. Henning. fe Mathis 



1992|). 



The first laboratory spectra used in astronomical silicate studies were of crystalline 
silic ates and natural gl asses (obsidian and basaltic glass, which contain some microlites; 



Pollack et al.lll973l ). Various st udies produced "amorphous" samp les through chemi 



e.g. 

cal vapor deposition fe. g.. iDavl Il979h . smokes (e.g.. iNuth fe Donnlll982l ). ion-irradiation of 



crystalline samples (e .g., Kratschmer fe Huffman 19791 ), l aser ablation of crystalline samples 



Dorschner et al.lll995l ; iJager et al. 



JScott fe Dulev!ll996h . and quenching melts to glass [e.g. 
20031 ). However, these techniques yield different results. For instance the peak positions, full 
width at half maxima (FHWMa) and ratios of the strengths of the 10 and 18 /im features 
vary between datasets even though the materials investigated are ostenibly the same com- 
position. Spectra from samples with the same reported composition vary, which should not 
occur if the samples had the same structures, and if all spectra were obtained under optically 
thin conditions. A detailed comparison between existing and new laboratory data are given 
in § [3j The sample preparation techniques vary widely and lead to a range of disordered 



2 Usually these are called optical "constants", but they are wavelength-dependent quantities, so we prefer 
•'functions" 
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structures. Unfortunately we do not have sufficient information on the physical structure or 
chemical characteristics of previously studied samples to determine the effects quantitively, 
but a qualitative analysis presented here highlights the need to make such sample information 
available. 

The data used m ost commonly for mo deling; astronomica l environments are synthetic op- 
tical fu nctions such as lDraine fc Led ( 119841 hereafter DL) and lOssenkopf. Henning. fc Mathis 
(119921 . hereafter OHM). These data are favored over laboratory spectra because they have 
broad wavelength coverage, which is not true of individual laboratory datasets. However, 
these synthetic functions were produced using compilations of laboratory spectral data and 
astronomical observed dust opacities from which new optical functions were calculated. The 
derived optical functions were then modified specifically to match the astronomical observa- 
tions. For instance, in the iDLl data, the 9.7 fim feature is entirely derived from observational 
opacities, while the NIR-NUV section of their optical function uses laboratory data from 
crystalline olivine studies, and the FU V/X-ray region uses laboratory data for crystalline 
alumina (AI2O3). Both lDLl and lOHMI blend laboratory and astronomical data and their op- 
tical functions will match some spectral features, and can be used for comparison of optical 
depths between different dusty environments. However, they do not represent real solids and 
thus cannot be used to determine the true nature of dust in space, how it varies spatially or 
temporally, or why. 



1.3. A brief guide to the structure of silicate minerals, glasses and the glass 

transition 

The basic building block of silicates is the SiO^ - tetrahedron. These can be linked 
in a framework, with each oxygen shared between two tetrahedra (e.g. SiC>2 minerals and 
feldspars), or they can be linked in chains (e.g. pyroxenes such as diopside [Di; CaMgSiaOe], 
enstatite [En;MgSiOs]) or they can be isolated tetrahedra (e.g. the olivines series: forsterite 
[fo; Mg 2 Si0 4 ] to fayalite [fa; Fe 2 Si0 4 ]). In all cases, the non-shared oxygens (known as non- 
bridging oxygens or NBOs) are charge-balanced by other cations (e.g. Mg 2+ , Fe 2+ , Ca 2+ , 
Na + , etc). The Ca-Mg-Al silicates that are expected to form in circumstellar environments 
are dominantly of pyroxene and olivine composition. 

Within each mineral group, solid solution allows compositions to vary between end- 
members. In the case of similarly-sized cations the solution may be continuous e.g., enstatite 
to ferrosilite (fs; FeSiOs) and the olivine series. Given the availability of other cations, end- 
member pyroxenes rarely occur in terrestrial or meteoritic samples. Where the substituting 
cation is a different size, solid solution is more limited, for example between enstatite and 
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diopside. In such cases, solid solution becomes less extensive at lower temperatures, so that 
cooling may lead to "exsolution" of two different pyroxene compositions from an initially 
homogeneous high-temperature solid. 

Minerals occur as crystals that possess long-range order, with very narrow distributions 
of bond angles and lengths. This leads both to anisotropy (properties varying with crystal 
orientation) and narrow spectral features. Silicate glasses are the "frozen" structural equiv- 
alents of liquids, possessing short-range order (so that local charge-balance is conserved for 
example) but lacking the long-range order that gives rise to symmetry and anisotropy in 
crystals. Glasses are therefore both isotropic and have broad spectral features. 

The basic structural unit of silicate glasses and melts is the SiO^ - tetrahedron, as is 
the case in crystalline silicate minerals. Oxygens linking tetrahedra are known as bridging 
oxygens (BO), while non-bridging oxygens (NBO) are coordinated by metal cations, which 
are termed network-modifiers in this role. Tetrahedral cations (T) include not only Si 4+ , 
but also trivalent cations such as Al 3+ and Fe 3+ ; these must be charge-balanced by other 
cations (usually alkalis or alkaline earths; Fig. [2]). The degree of polymerization of a melt or 
glass can be summarized by the ratio of NBO/T, which can range from (fully polymerized, 
e.g. Si02) to 4 (fully depolymerized, e.g. Mg 2 Si04). In general, more polymerized melts 
are more viscous and have higher glass transition temperatures. On quenching a melt, 
its structure is "frozen in" at the glass transition if cooling is rapid enough to prevent 
crystallization. The glass transition is actually an interval, often approximated by the glass 
transition temperature (T g ) which is usually taken to be the temperature at which the 
viscosity is 10 12 Pa.s (T12). Rapid cooling from a given temperature preserves the network 
present in the liquid at that particular temperature. Because of this behavior, glasses of 
any given composition can have subtle differences in structure that depend on cooling rates. 
The temperature at w hich the glass ha s the s ame structure as the melt is called the Active 



temperature (Tf). See lMysen fc Richetl (120051 ) for a comprehensive review of melt structure 
and properties. 

Two issues whose importance will be discussed in the current work are the role of water 
and the oxidation state of iron. At low water contents (less than about 1 wt. % total H^O) , 



water dissolves in silicate glasses almost exclusively as hydroxyl (OH - ) ions (IStolperlll982l ). 
and acts as a network modifier (Fig. [2]). Compared to other modifier oxides such as Na 2 
and MgO, wat er has a more drama tic effect in reducing melt viscosity and glass transition 



temperatures (IDingwell et al.l 119961 ) . Iron can play the role of network modifier (octahedral 
Fe 2+ or Fe 3+ ) or network-forming cation (tetrahedral Fe 3+ ). Consequently, the oxidation 
state of an iron-bearing glass or melt has a significant effect on its structure and properties. 



From the perspective of cosmic dust formation, the glass transition temperature is es- 
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sentially the temperature above which a given composition should form as or convert to 
crystalline solids, w hereas solids formed below this temp erature will be glassy if they cool 
sufficiently rapidly (IRichet et al.l Il993l ; ISpeck et al.l 120081 ). More depolymerized melts are 
more difficult to quench, and melts less polymerized than pyroxenes (NBO/T > 2) typi- 
cally require extreme quench rates (100s Ks" 1 ) using methods such as containerless laser 
processing to achieve truly glassy samples (ITangeman et al.ll200ll ). For a given composition, 
faster cooling rates result in a higher T g . This dependence can be determined by differential 
scanning calorimetry using different heating and cooling rates, and is used to determine the 
cooling rate of natural lava samples (jWilding et al.lll995l ). 



For depolymerized silicates (e.g. olivines and pyroxenes), if glassy grains form they must 
do so below their T g because the cooling rate required for quenching to a glass is extremely 
rapid. Highly polymerized silicat es (e.g. silica, obsidian) can be cooled more slowly, over 
hours or days, and not c rystallize (ISchairer fc Bowenlll956T) . However, the cooling timescales 
(months) determined by lSpeck et al.l (120081 ) for AGB star circumstellar shells are sufficiently 



long as to preclude the preservation of glassy/chaotic solids that form above T g , because 
annealing timescales are shorter than those for cooling for all but the most polymerized 
silicates. 



1.4. The need for new data 

Modeling of silicate dust in space has been limited by t he available laboratory data 



The influence of various model parameters was investigated by I Jones fc Merrilll (119761 ) . who 
found that using so-called "clean" (i.e. pure magnesium) silicate grains to model the observed 
9.7 //m features did not yield a good fit due to the lack of absorption by these grains in the 
visible and near-IR. They also found that just mixing in more absorbing grains did not solve 
the problem. This led to the suggestion that the grains responsible for the 9.7 //m feature 
are "dirty" silicates, i.e. Mg-silicates with impurities introduced into the matrix giving more 
opacity in the optical and near-IR. 

It is known that NBO/T (polymerization) affects the spectra of amorphous silica tes suc h 



that the peak position of the 10 /im feature shifts redwards as NBO/T increases (e.g. lOHMl ). 



Aluminium (Al) is a n etwork former and consequently Al content strongly affects NBO/T. 



Mutschke et al.l (119981 ) suggested that Al may be an important component of silicates in 
space that could explain why previous laboratory spectroscopic studies failed to match ob- 
servational data. Other cations may be equally important. Ca 2+ and Fe 2+ both substitute 
for Mg, while Fe 3+ will substitute for tetrahedral site (e.g. Si 4 + or Al 3+ ). Therefore, the 
oxidation state is another important variable in addition to elemental substitutions. 
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Existing laboratory spectral data for "amorphous" silicates were produced using samples 
that are not sufficiently well-characterized to allow astronomers to interpret their observa- 
tions without ambiguity. Here we present new laboratory spectra of several silicate glasses 
of astronomical relevance, and discuss compositional factors that influence their spectral fea- 
tures. We compare these new data with those previously available for "amorphous" silicates 
and discuss how these samples compare to successfully-applied synthetic optical functions. 
We find that the synthetic spectra cannot be well matched by the conventionally assumed 
glassy silicate composition and discuss whether astrophysical silicates need to be truly glassy. 



2. Experimental Methods 



2.1. Sample selection 

The bulk composition of silicate Stardust lies somewhere between pyroxene (M2S12O6) 
and olivine (M^SiO^), where M indicates metal cations, with Mg and Fe being the most 
abundant. The cosmi c Mg/Si ratio is predicted to be ~1.02, while Fe/Si ~0.84 (e.g., 



Lodders fe Fegleyl Il999l ) . However, most Fe is ex pected to combine with S, Ni, Cr, Mn and 
other siderophile elements into metallic grains ( jGail fe Sedlmayrl Il999l ; lLodders fe Fegley 
19991 ). This partitioning is reflected in the Earth, where most iron resides in the metal- 
lic core while the mantle is dominated by magnesium-rich silicate with Mg/Fe ~9. The 
predicted bulk cosmic silicate would then be close to MgSiOs, with minor amounts of iron 
leading to an atomic (Mg+Fe)/Si ratio slightly greater than 1. Determining the spectra of 
various olivine and pyroxene glasses is therefore of critical importance for identifying the 
silicate mineralogy. The focus of this study is glass compositions for which data already 
exist in the astronomy literature i.e. predominantly olivines and pyroxenes (See Tableland 
references therein). Mg-rich endmembers are forsterite and enstatite, respectively. Another 
mineral that is commonly discussed in astromineralogy is diopside, which is als o a pyroxene. 



Diop s ide has been invoked to explain observed crystalline silicate fe atures (e.g. iDemyk et al. 



2000l ; iKemper et al.ll2002t iHony et al.ll2009l ; lOnaka fe Ok ada 2003|) and appears in the clas- 



Tielensl ll99oT see Fig [3]). Furth e rmore 



sic condensation sequence for dust formation (e.g. 
alu minous diopsid e form ed in the experimental condensation study by lToppani et al.l ( 120061 ). 
and IDemyk et al.l (120041 ) showed that crystalline diopside grains can easily be amorphized 
by heavy ion irradiation. 

To complement the olivines and pyroxenes (Table [1]) we include four samples from the 
melilite series, whose endmembers are gehlenite (Ge; Ca2Al 2 Si07) and akermanite (Ak; 
Ca 2 MgSi 2 07) . gehlenite is predicted to be among the first-formed silicate grains (e.g., 
Tielend Il990l ; lLodders fe Fegleyl Il999l . see Figure [3]), and the major repository for calcium 
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and aluminum in dust, whereas pyroxenes are predi cted to be among the most abundant 
grains, and the major repository for Mg and Si (e.g 



Tielens 


1990; 


Lodders & Feglev 


1999 



Gail fe Sedlmayrlll999l ). Furthermore aluminium a nd calcium are highly depleted from the 



gas phase and are assumed to be included in dust (jWhittetlll992l ). Aluminum-rich silicates 



like gehlenite and other melilites are major constituents of Calc ium- Aluminum- ric h Inclusions 



(CAIs) and gehlenite has been identi fied in red super giants ( iSpeck et al.ll2000|) and Active 



Galactic Nuclei (AG N) environments (IJaffe et al. l2004 and in meteorites (IStroud et al.l 2008: 
Vollmer et al.l 120071 1 . Furthermore, a re cent Type lb super nova (SN2005E) has been shown 

, Inde ed almost 50% of its ejecta 
( 120071 ) reported laboratory spec- 



to contain more calcium than expected flPerets et al.l 2010 
mass (or >0.1M Q ) is attributed to calcium. IChihara et al 



tra of crys talline melilites e very 10% along the solid-solution join between akermanite and 
gehlenite. iMutschke et al.l (119981 ) reported laboratory spectra of two glasses: end-member 
gehlenite, and ak50-ge50. Given that the cosmic Mg/Ca ratio is ~16, it seems likely that 
melilites could have substantial akermanite contents. The difference in the structure of the 
glasses is profound: NBO/T is 0.67 for gehlenite and 3.0 for akermanite, while Al/Si is 2 for 
gehlenite and for akermanite. 



Both theoretical models (ILi fe Lill2009l ) a nd observations (IJones fe Merrillll 19761 ) suggest 
that some iron is incorporated into silicates. iJager et al.l ( 119941 ) presented spectra for two 
silicate glasses containing the most abundant dust forming elements (i.e. Mg, Si, Fe, Ca, Al, 
Na). However, iron in their sample was partially oxidized (FeO/Fe203 ~1) which leads to 
problems in interpreting the spectrum (see § ll.3p . With this in mind we present a sample 
we call "Basalt" (Nao.09Mgo.62Cao.69Feo.39TioxoAlo.06Si2.i6O6)) which contains ferrous iron 
(Fe 2 +). 

While iron has been invoked to explain opacity problems, most iron is expected to 
combine with other siderophile elements to form meta l or metal sulfide grains rather than 
silicate ( jGail fe Sedlmayrlll999l : lLodders fe Fe gleyl 1 1 99 91 ) . Consequently, we have synthesized 
an iron-free silicate glass using cosmic abundance ratios for Mg, Si, Al, Ca, and Na, yielding 
(Nao.ioCao.i2Mgi. 86 )(Alo.i8Sii.84)06. This sample does not include volatile elements or iron 
and is quite close to en statite in composition. This "cosmic silicate" is designed to test 
the IStencel et al.l (Il990h hypothesis that dust forms as chaotic solids with the elemental 
abundances of the gas. 



Pollack et al.l (119731 ) provided spectra of obsidian, a naturally occurring glassy silicate. 
Consequently we include obsidian glass in our sample, in part because this is the origin of 
the attribution of the ~10 /zm feature to amorphous silicates. It is also useful for studying 
the effect of silicate structure and composition on spectral parameters because it is different 
from the more commonly assumed olivines and pyroxenes. 
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Finally, we include silica (Si0 2 ) glass in our sample. Like obsidian, the structure of 
Si0 2 is significantly different from olivines and pyroxenes and thus provides potential insight 
into the effects of structure on silicate spectral features. Furthermore, silica dust grains 
have been invoked to explain observed astronomical spectral features in both evolved stars 
(ISpeck et al.ll2000l ) and in young stellar objects (e.g. ISargent et alll2006l ). 

The samples investigated here are listed in Table |2j 



2.2. Sample synthesis and preparation 

Samples designated with "synthetic" in their name are generated from mixtures of 
reagent-grade oxides and carbonates, providing glasses with compositions like those of end- 
member minerals. The "cosmic" silicate was also synthesized in this way. In contrast, 
samples designated as "remelt" were generated by melting natural mineral samples. Conse- 
quently the "remelt" -samples have compositions whose additional components reflect the im- 
purities^ found in natural crystalline samples of the relevant minerals. Using both synthetic 
and remelted samples helps to demonstrate the effect and importance of small compositional 
variations in silicates. 

Synthesis of silicate glasses from oxide and carbonate starting materials, was undertaken 
at the MU experimental pet rology facility and the procedures are described in detail by 



Getson fc Whittingtonl (120071 ). Melilite glasses were prepared by fusion in Pt crucibles and 
quenched by pouring into graphite molds (slow cooling), or on to a copper plate (faster 
cooling for less polymerized compositions). Glasses in the pyroxene series (including "cosmic" 
compositions) vary in their quenchability; diopside (CaMgSi 2 6 ) is an excellent glass-former, 
while enstatite (MgSiOs) crystallizes extremely rapidly. 

Glass of forsterite (Fo, Mg 2 Si0 4 ) or olivine (Mg 2a ,Fe 2 _ 2 .j;Si04) can only be formed in 
the laboratory under special conditio ns. Specifically, Fo glas s in the form of 50 to 200 /zm 



diameter beads was first produced by iTangeman et al.l (120011 ) at Containerless Research by 



suspending small particles in argon gas and melting/quenching by pulsing with a laser. In 
addition to providing rapid cooling of ~700 K/sec, the lack of container promotes crystal-free 
glass formation. Commercially prepared samples of Fo were purchased from Containerless 
Research, Inc. Iron-bearing olivine samples were not available. We were able to prepare 



3 Unlike minerals, glasses do not have to have well-defined formulae. Some studies describe such non- 
stoichiometric glasses as having large quantities of impurities, but these non-mineral-end-member composi- 
tions are simply what the glass is made of. 
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"basalt glass" at Washington U. by flash melting rock on a glassy carbon substrate in a 
vacuum chamber with a CO2 laser followed by rapid cooling. When applied to iron-bearing 
olivine or fayalite (Fe 2 Si0 4 ) crystals this approach failed to produce glass. 



2.3. Chemical Analyses for Sample Composition Determination 



Samples were characterized by wavelength dispersive analysis (WDS) using standard 
procedures on the JEOL-733 and JXA-8200 electron microprobes at Washington University 
using "Probe for Windows" for data reduction (see http://www.prob esoftware.com/). The 
measured data were corrected using CITZAF after lArmstrona (119951 ). Oxide and silicate 
standards were used for calibration (e.g., Amelia albite for Na, Si; microcline for K; Gates 
wollastonite for Ca; Alaska Anorthite for Al; synthetic fayalite for Fe; synthetic forsterite for 
Mg; synthetic TiO"2 for Ti; synthetic Mn-olivine for Mn; synthetic C^Os for Cr). Microprobe 
analyses are given in Table [3J and the resulting compositions are given as chemical formulae 
in Table [5J Table [3 ] also lists water c onten ts, which we determined from near-IR spectra 
using the method in iHofmeister et al.l (120091 ) . 



2.4. Viscosimetric Determination of Glass Transition Temperature 

As discussed in § ll.3[ the glass transition temperature (T g ) is essentially the temperature 
above which a given composition should form as or convert to crystalline solids, while solids 
formed below this temperature should be glassy. We determine T 12 values as a proxy for T g for 
our silicate samples to provide an upper limit on temperature for models of glass formation in 
space. T g depends on composition and cooling rate. Viscosity was measured over a range of 
temperatures using a T heta Instruments Rheotro nic III parallel plate viscometer, following 



procedures described by lWhittington et al.l ( 120091 ). Viscosity is calculated from the measured 



longitudinal strain rate, known load and calculated instantaneous surface area, assuming 
perfect slip between sample and plates. The measurements are interpolated to find T12 with 
an uncertainty of less than 2 K. 



2.5. Spectroscopy 

Room temperature (18-19°C) IR absorption spectra were acquired using an evacuated 
Bomem DA 3.02 Fourier transform spectrometer (FTIR) at lcm -1 resolution The accuracy 
of the instrument is ~0.01 cm -1 . Far-IR data [y < 650cm _1 ) were collected for five samples 
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using a SiC globar source, a liquid helium cooled bolometer, and a coated mylar beam- 
splitter. Mid-IR data {y = 450-4000 cm -1 ) were collected for all samples using a SiC globar 
source, a liquid nitrogen cooled HgCdTe detector, and a KBr beam-splitter. The spectra were 
collected from powdered samples pressed in a diamond anvil cell (DAC). The empty DAC 
was used as the reference spectrum which allows ref lections to be removed. The methodology 
is described by iHofmeister. Keppel fc Speckl ( 120031 ) . 

Interference fringes exist in many of our spectra because the diamond faces are parallel 
and are separated by a distance within the wavelength range studied. The spacing is larger 
than film thickness. Fringes are associated with a sideburst in the interferogram. Due to 
mathematical properties of Fourier transforms, the fringes are convolved with the spectrum 
and therefore do not affect peak shape or the parameters used to describe the peak (position, 
width and height). 



3. Comparison of new laboratory spectra with previous laboratory and 

synthetic spectra 

3.1. New laboratory spectra 

The new laboratory spectra, shown in Figure HI are the highest resolution spectra of sili- 



cate glasses to date. These data are available online from http:/ /galena. wustl.edu/~dustspec/idals. html 
The main observable parameters of astronomical spectra are the peak position of the ab- 
sorption/emission feature s at ~ 10/mi and ~ 18/im; their FWHMa and the ratio of their 



strengths (see e.g., lOHMI and references therein). For instance, the ratio of band strengths 
between the 10 and 18 /im features in observations varies markedly, as do the peak positions. 
Therefore we have extracted the most important spectral parameters from our data using 
the NOAO onedspec package within the Image Reduction and Analysis Facility (IRAF). 
The spectral parameters (the peak position, barycentric position, full width half maximum, 
FWHM, and equivalent width) were determined for the ~ 10 and ~ 18/xm features. Figure H] 
shows the barycentric position of the ~ 10/im feature for all samples, while Figure |5] shows 
the barycentric positions of both the ~ 10 and ~ 18/mi features for the five compositions for 
which far-IR data were collected. The barycentric positions, along with the peak positions 
and FWHMa from absorbance measurements are listed in Table HI 

There are multiple terms and symbols for the various ways in which absorption of light 
by solids is described. In order to prevent confusion, we will explain precisely how each 
term we use is defined. This is particularly important for applying laboratory spectra to 
astrophysical studies because our spectra are initially in the form of absorbance (a), but we 
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typically use either absorption efficiency (Q- values), optical depth or extinction in astronomy. 



Transmittance, T (which is also referred to as transmissivity in iFoxl 120021 ) is defined as 
the ratio of the intensity of transmitted (/trans) and incident (I ) light, i.e. Itrans/^o- Similarly 
reflectivity, R and absorptivity, A are the ratios of absorbed to incident light and reflected 



to incident light, respective! 



ptivity 
iveljfl 



-^abs -^trans I) 



rcfi 



-^trans ^ -^abs 



-frefl 



T = 1-R-A 



For simplicity we will assume the reflectivity is negligible or has been accounted for (see 
Hofmeister et al.ll2009l . for how we can account for reflections) . Then, 



T = 1 — A 



When light passes through a solid the absorption can be expressed as: 



he 



-aL 



(2) 



where L is pathlength or thickness of the solid sample and a is usually called the absorp- 
tion coefficient, but is sometimes called opacity. In addition a = up, where k is the mass 
absorption coefficient, and p is the mass density. 



T 



-aL 



Absorbance, a, is the exponent in the decay of light due to absorption: a 
Optical depth t\ of an absorbing material is defined by: 



aL 



(3) 
npL. 



note that absorb ance is the reciprocal of the log of the transmittance, in contrast to the absorptivity 
defined here; see also lSpeck. Thompson fc Hofmeisterl (|2005j ) 
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I e~^ (4) 



(From, e.g. lGlassHl999f ). 



From equations [21 [3] and H] we see that the absorbance a is similar to optical depth T\. 
In order to compare astronomical data in which we have a wavelength-dependent optical 
depth we can use absorbance (which is how we present our transmission spectra in Fig. H] 
and Fig. [5] after accounting for surface reflections.) Now we can relate the absorbance (a) 
and absorptivity (A) via equation [3j 



A = 1 - e 



-a 



To compare with some astronomical observations we still need to extract a version of the 
laboratory data that is comparable to the absorption efficiency, Q-values. To get this we 
need to consider how Q-values are defined. For a non-blackbody dust grain we define an 
absorption cross-section C a b s as the effective geometrical cross-section of the particle once 
we account for it not being a blackbody: 

Cabs = Qabs x Y 

where T is the geometrical cross-sectional area of a dust grain. Now if we consider how the 
absorption cross section gives rise to absorption we get: 

^ = A = C, hs nL (5) 

where n is the number density of absorbing particles and L is the pathlength or thickness of 
the absorbing zone. 

n = T7—^ ( 6 ) 

M mo i x m H 

where M mo i is the molar mass of the solid and m# is the mass of a hydrogen atom. Combining 
equations [5] and [6] we get: 

_ A x M mol x m H m 

Vabs Y pL 
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It is clear from equation [7] that Q- values oc absorptivity such that: 

Qabs = (A 

where, 

_ M mol x m H 
^ ~ TpL 

Consequently, while the shape, peak position and FWHM of spectral features shown in Q- 
values will be identical to those for A, the absolute values depend on the pathlength and on 
the cross-section areas of a given grain distribution. The uncertainty in the thickness of our 
samples and in potential grain size makes the absolute Q-values uncertain. Consequently we 
normalize our Q- value spectra to peak at unity. These Q- value spectra are stacked in Fig [6] 
to demonstrate the shift in both barycentric position and FWHM with composition. Table H] 
also includes the FWHM in Q- value. 



3.2. How to compare disparate data sets 



Here we compare the new laboratory spectra with those previously published to distin- 
guish which factors are most important in determining spectral feature parameters. Previ- 
ously published data are available as complex refractive indices (n and k). To make a fair 
comparison between the many available datasets we converted our absorbance data to the 
wavelength-dependent imaginary part of the complex refractive index (k) for each sample. 
This absorption index k is chosen as the best comparison of different datasets, as it does 
not depend on grain size s or shapes and can be extracted equally well from transmission or 
reflectance data (see e.g.. lHofmeister et al.ll2009l ; lFoxll2002[ ). Comparing absorption efficiency 
Q a bs requires assumptions about grain shape, which have been shown to affect the shapes 
and positions of spectral features (e.g. iMin et al.ll2003l ; iDePew et al.l 120061 ) . The conversion 
of the absorbance spectrum to fc-values uses: 



k = 2.303a/ (d4%v) 



where a is the absorbance, d is the sample thickness (in cm) and v is the frequency (wavenum- 
ber) in cm -1 . There is some uncertain ty in the measurement of th e sample thickness which 
is estimated to be 0.5 - 1.5 //m thick ( IHofmeister &: Boweyl 120051 ) . Consequently, for com- 
parison our data has been normalized to peak at a k- value of 1. In Figure [7] we compare 
spectra of forsterite, enstatite, gehlenite and akermanite composition "amorphous" silicates. 
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3.3. Comparison of new glass spectra with existing laboratory data. 



Various synthesis techniques are associated with samples studied in the laboratory, as 
described in § 11.21 and listed in Table [H Considering, for example, forsterite composition 
samples shown in the upper left panel of Figure Ul it is clear that the spectral parame- 
ter s vary even for ostensib l y the same composition. The ion irradiation technique used 
by iKratschmer fc Huffman! (119791 ) apparently does not produce a fully amorphized sample, 
as this spectrum is closer to that of crystalline forsterite, which was their st arting material. 
The sa mple from lDayl (119791 ) was produced via chemical vapor deposition; the lScott fc Duley 
(119961 ) samp le was produced by laser ablation of crystalline samples; and the sol-gel method 
was used by IJager et al.l (120031 ) . 



The spectra of three different samples produced by chemical vapor deposition, laser 
ablation and sol-gel techniques are similar. However, although these samples may be "amor- 
phous" they are not necessarily glassy. The samples generated by these three techniques may 
represent chaotic silicates rather than the glassy silicates investigated herein. The difference 
between previous samples and those presented here is most likely a combination of density 
and porosity. Glasses should be less porous and denser than chaotic silicates. 



Among the enstatite composition samples, those from iDayl (Il979l ) and IScott fc Duley 



(11996) are, again, sim i lar to each other. Samples produced by melting and quenching (e.g. 



ours, iDorschner et al.l Il995l ) show some variability in the spectral features. However, the 
peak positions and FWHMa are similar (see Tables [J and [2]) . The difference in breadth may 
result from different cooling rates and Active temperatures. IJager et al.l (120031 ) investigated 
amorphous MgO-Si02 solids prepared using the sol-gel method. The spectra vary markedly 
with changes in SiO"2 content, and hence with polymerization (Fig. [H]). However, these 
samples also contained a reported 0.8 to 1.2 wt.% H 2 0, which equates to ~3 mol.% H 2 0, all 
dissolved as network- modifying hydroxyl (OH~) ions fse e §ll.3|). Glasses quench ed from melts 
at atmospheric pressure (e.g. those presented here and in lDorschner et al.lll995l ) contain much 
less water, typically 0.02-0. lwt.%. In previous laboratory studies it has been suggested that 
the changes in peak positions are due to differing silica contents, and hence polymerization 
states (NBO/T), i.e. forsterit e has a redd er peak position than enstatite, which in turn 



is redder than silica (see e.g., lOHM] Il992l . and reference therein). This is demonstrated 



in Fig. [HI Thus, even the modest levels of water remaining in samples prepared using the 
sol-gel method will affect the structure of the silicate and thus its spectrum. This makes 
precise interpretation of spectra from sol-gel samples difficult, especially of the more silica- 
rich compositions whose structure will be the most affected by the incorporation water. 
Therefore the sample / material structur e is not well known. Spectral dif ferences between 



the sol-gel MgSi0 3 of IJager et al.l (120031 ) and the quenched MgSi0 3 glass of IDorschner et al. 
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( 119951 ) emphasize the point, especially around the 10/im feature (Fig. [7]). 



The bottom left panel of Fig. [7J compares spectra for the gehlenite, the Al-rich endmem- 
ber of the melilite series. All data were pro duced by the mel t -quen ch method. Gehlenite from 
this study closely matches spectra of the iMutschke et all (119981 ) sample, while the sample 
we generated by melting a natural crystal of gehlenite (designated gehlenite remelt) shows 
significant differences which can be attributed to the deviation from stoichiometry which 
results in large shifts in both Al/ (Al+Si) and NBO/T (see Table [21) • The bo t tom r ight panel 
compares the mid-composition melilite (Ak50Ge50) from IMutschke et al. ( 1998 ) with our 
end-member synthetic melilites. The peak position of the IMutschke et al. I dl998h melilite is 
similar to our Al-rich endmember, while the feature shape is closer to that of our Mg-rich 
endmember. 



Mutschke et all (119981 ) discussed the low contrast feature in the 12-16 /im range in 
various aluminous silicates. This feature moves from ~ 14.5 /im in gehlenite towards shorter 



wavele ngths, with the akermanite feature peaking closer to 13.5 /im (Fig. [7]). IMutschke et al. 



(119981 ) did not investigate the akermanite-rich melilites because their focus was on the effect 
of aluminium, but their findings , and those seen here may pertain to the carrier of the 



observed "~ 13 /mi" feature (e.g. ISloan et al.ll2003l . and references therein). 



Figure [9] shows a comparison bet ween our "cosmic silicate" , "basalt" and the "dirty 
silicate" produced by lJager et all ( 11994J ). The 10 /im features are very similar in peak position 
and FWHM, but the ratio of the ~10 and ~ 18 /im features vary. This comparison shows 
that samples containing several elements give rise to very similar 10 /im spectral features 
even though they differ in Fe/(Mg+Fe), oxidation state, NBO/T and other compositional 
parameters. 



I n single crystal silicates, the Mg/Fe ratio affect s the s pectral features (e.g.. iKoike et al. 



20031 ; iHofmeister fc Pitman! 120071 ) . iDorschner et all (119951 ) investigated the effect of Mg/Fe 
ratio on the spectra of (Mg,Fe)SiC>3 glasses produced by quenching melts in air. The peak 
positions shift slightly, and the ratio of the 10 /im and 18 /im feature heights varies markedly 
(Fig. |Sj). Substitution of Ca for Mg produces substantial broadening of the 10 /tm feature, 
consequently the peak (and barycenter) shift redward. The viscosity and glass transition 
data in Table [2] show that small impurities can have important effects on melt properties 
(and structure). For example, compare synthetic and remelted akermanite and gehlenites. 
This is why it is important to consider Ca in pyroxenes, in addition to the En-Fs series. 
Indeed based on the spectra presented here, the effect of calcium substitution is larger than 
that of iron. 



Dorschner et all (119951 ) reported that the FeO/Fe203 ratios of their samples were 
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Whether molar or weight ratio (unspecified in their paper), these samples do not have py- 
roxene stoichiometry, and thus have a structure differing from that of our pyroxene glasses. 
Using their reported compositions, and assuming all Fe 3+ acts as a network former, calculated 
NBO/T values ranges from 2.1 for En95 to ~0.8 for En50 and En40 glasses, considerably 
lower than the value of 2.0 for a true (Mg,Fe)SiC>3 composition. Furthermore, the effect 
of tetrahedral Fe 3+ on neighboring Si-0 bonds, which give rise to the 10 /im feature, may 
be sufficient to produce marked changes in the shape and peak position of the feature even 
if Mg-Fe substitution does not. Therefore existing lab data on glasses do not allow astro- 
nomical spectra to be interpreted reliably in terms of either M 2+ /Si or Mg/Fe ratio, which 
are important tools for discriminating between competing dust formation models. We will 
address the roles of oxidation state and Al-content in silicates in future papers. 



3.4. Comparison with Synthetic Spectra 



As discussed in § 11.21 the most popular silicate spectral data used in astronomy are the 
synt hetic o ptical functions (complex refractive indices/complex dielectric functions) from lDL 
and lOHMI . In the spectral region s con sidered here, both groups derived di electric functions 
from astronomical observations. iDLl used observations of the ISM while lOHMI produced 



two sets of optical constants designated as warm O-deficient and co ld O-rich. The warm 
O-deficient is intended to match circumstellar dust features where the iNuth &: Hechtl (119901 ) 
predicted that non-stoichiometric silicates form in the outflows. The cold O-rich silicate 
is meant to represent dust formed in molecular clouds, where the cool temperatures and 
slow dust formation lead to sto ichiometric composit ions. For the circumstellar (O-deficient 
warm), opacities are derived by lVolk fc Kwokl (119881 ) based on averaging approximately 500 
IRAS LRS spectra of evolved stars. However, this approach is less than ideal. It is well 
known that the observed "silicate" features vary markedly from source to source (see e.g. 
Nuth &: Hechtl 11990c iMutschke et all Il998t ISpeck et all |2000| ) and thus averaging can smear 
out such differences and give rise to an opacity that is close to matching many objects and 
actually matches none. 

To provide a fair comparison, para meters for these synthetic silicate spectra have been 



included in Table [TJ Whereas the IDLl spectra provide a reasonable match to the features 
observed in the ISM, the spectra from lOHMI provides a nar rower ~ 10/xm featu re which 



matches the circumstellar silicate featu res more closely (e.g. ISargent et all |2010| . and ref- 



erences therein). However, neither IDLl nor lOHMl provide perfect matches to the observed 



astronomical mid-IR features; nor do they reflect the diversity of features seen in astronom- 
ical environments. Although synthetic spectra do not match astronomical data perfectly, 
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they are widely used and come close to matching many astronomical observations. Conse- 
quently, comparing our new laboratory data to these synthetic spectra allows us to assess 
whether any true glasses show promise as carriers of observed silicate features. 

Figures [TOltofT^lshow how the spec tral proper ties of the new laboratory samples compare 
to those of synthetic silicates from Idl] and OHM . In all cases the comparison is made with 
the absorption index k, which is the imaginary part of the complex refractive index, and was 
chosen to avoid complications arising from grain shape and grain effects (see § I3.2p . None 
of the samples presented here precisely match the sp ectr al features of the synthetic silicates. 
In particular, the 10/zm feature in the spectra from iDLl is too broad to be matched by any 
of the laboratory samples, with the possible exception of gehlenite (both "synthetic" and 
"remelt"; see Fig. [TT] bottom row). However, the relative cosmic abundances calcium and 
aluminium make this an unlikely attribution for the ISM dust. It is usually assumed that 
astronomical silicates are largely comprised of "amorphous" Mg-rich olivine and pyroxene 
compositions. Therefore, we compare the interstellar feature as represented by lDLl with such 
compositions in Fig. [T2TIT31 The spectra of the four Mg-rich pyroxene-like glasses (Synthetic 
Enstatite, Enstatite Remelt, Cosmic Silicate and Basalt) h ave g rossly similar spectra. In all 
four cases, the laboratory data match the blue side of the IDLl feature, but are too narrow 
and fail to match the red side. This is also true for the broader diopside feature (Fig. [T2|) . 
While forsterite composition glass has a red der f eature than the pyroxenes, it is still too 
narrow and still not red enough to match the IDLl 10/zm feature. 



If glasses do produce the astronomical 10/zm band, the excess breadth n eeds accounting 
for and may suggest processing of dust in the ISM (INuth fc Hechtl Il990l ). The physical 
changes to the dust that give rise to the broader interstellar feature remain unknown, but 
there are numerous interpretations. In circumstellar environments the red-side broadening 



of th e 10 /im feature has been interpr eted as being due to oxide inclusions (e.g..lSpeck et al. 



in grain shape (e.g.. iMin et al 



20001 ) ; increasing cr ystallinity (e.g., ISylvester et al.l Il998t iBouwman et al.l 120011); changes 



20071 ); or grain porosity (e.g., iHenning &; Stognienkol Il993 



Voshchinnikov fe Henningj 120081 ). The silicate feature is fairly constant and broad for the 



diffus e ISM, but varies more and is narrower for molecular clouds (see e.g. Ivan Breeman et al. 



Chiar & Tielens 


2006; 


Chiar et al. 


2007 



McClurd 120091 ; Ivan Breeman et al.ll201ll ). but the results of those studies depend on the 



optical properties input into their models. 



The spectra of O HMl have been used somewhat successfully in modeling circumstellar 
silicate features (e.g. ISargent et al.l |2010| . and refer ences t herein). Of the new laboratory 
glasses, only forsterite comes close to matching the lOHMl synthetic spectra (Fig. [Hj), but 
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there is still a problem with the breadth and redness of the feature. Most other compositions 
have 10/zm features that are too blue. However, 10 /im silicate emission features seen in AGB 
stars spectra and o bserved by ISO and IRAS are consistently sl ightly bluer than those in the 
synthetic spectra (jwheeler fc Speck 2007 ; Wheeler et al 2011 ). Consequently our samples 
may provide better matches than the OHM synthetic spectra. 



4. Discussion 

4.1. Amorphousness at different scales: Glasses vs. Nanocrystalline solids 



We have shown that differences exist between spectra of amorphous samples of osten- 
sibly the same composition due to structural differences arising from synthesis methods. 
In this section we discuss the range of structures that may be considered amorphous, and 
other factors such as composite or polycrystalline grains that may produce different spectral 
features to those obtained from studies of single crystals. The range of structures is shown 
schematically in Fig. [TJ 

Disordered (but not glassy) silicates may form at temperatures above T g by ion bom- 
bardment of initially crystalline materials. Depending on the extent of damage done by ion 
radiation, an initially crystalline sample may be amor phized and could be indistinguishable 
spectroscopically from a truly glassy sample (see e.g. iDemyk et al.ll2004l ). As discussed in 
§ 11.21 limited heavy ion irradi ation may not completely amorphize a crystalline sample (e.g. 
Kratschmer fc Huffman! 1 19791 ). Consequently ion- irradiation of initially crystalline material 
should lead to a continuum of structures ranging from perfect crystals to completely disor- 
dered, with the accompanying range of spectral features. Furthermore IDemyk et al.l (120041 ) 
showed that ion irradiation leads to more porous samples than simple splat-quenched glasses. 
While lower densities should not directly affect the spectra, voids give rise to extra reflec- 
tions which increases extinction and feature widths through increased scattering even though 
absorption is unchanged. The structures of ion-irradiated crystals may be more similar to 
smokes than glasses because of ion damage (see Fig. [p. 

Another potential carrier of the observed silicate spectral features is polycrystalline sili- 
cate. Spectra of single composition (monomineralic) polycrystals will be affected by multiple 
sc attering whic h shou ld broaden and smear out the features as demonstrated for hematite 



in 



Marra et al.l (120111 ). During annealing and crystallization of an initially glassy grain the 



final structure would depend strongly on composition. For a glass of pure enstatite composi- 
tion, one may expect either a single crystal of enstatite or a mono-mineralic polycrystalline 
agglommeration (if crystallization starts at more than one point in the grain.) However, 
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since glassy grains are expected to form to include all the atoms in the outflowing gas, frac- 
tional crystallization is possible. For instance, our "cosmic" silicate has a bulk composition 
close to enstatite, but contains a significant amount of other elements. Consequently we 
might expect an annealed sample of this composition to be poly-mineralic polycrystals with 
a large abundance of enstatite crystals. Depending on the spectral contributions from other 
constituents of the polymineralic grain, we may not perceive sharp crystalline features in 
the spectrum. This is beyond the scope of the present work but will be investigated in the 
future. 

Given that there are multiple mechanisms for the formation of both amorphous and 
crystalline silicate grains, and that processing can lead from one to the other and vice versa, 
it is possible that many types of grains represented on Fig. [T] are found in space. 



4.2. Potential application to astronomy 



AGB stars present an interesting environment in which to study dust formation be- 
cause they are relatively benign and the stability of CO molecules simplifies the chemistry. 
Whereas many AGB stars exhibit the ~ 10/mi feature, this feature varies in p eak position 



and shape from object to object, and even t emporally w i thin a single object f lSpeck et al. 



20001 ; ISloan et al.112003 : Monnier et al.lll998l ). ISpeck et all (120001 ) s howed that s e veral Galac- 
tic objects show silicate features peaking as short as 9.2 /zm, while ISpeck et all ( 120081 ) found 
a very red 10 /im absorption feature in the spectrum of an obscured AGB star. While the 
synthetic spectra are commonly used to model the observed features, they rarely match the 
details of peak position and width. The new laboratory data presented here provide a frame- 
work to interpret t he observed v ariati ons in AGB star spectral features. For example, the 
red feature seen by ISpeck et all ( 120081 ) can be best matched by something forsteritic, while 
the observed 9.2/im feature needs very silica-rich dust. 



Recent studies bv IChiar fc Tielensl (I2006f ): IChiar et al.l (120071 ) . McClurd (120091 ) and 
van Breeman et all ( 1201 ll ) have shown that the shape and peak position of the classic 10 fim 
interstellar silicate absorption feature varies depending on the line of sight. Whereas the 
~ 10/mi feature remains the same for all diffuse lines of sight, its shape and position varies 
once the line of sight includes a molecular cloud. This has been attributed to a combination 
of grain agglomeration and ice mant le fo rmation within the molecular clouds (op. Cit.). The 
diffuse ISM is well characterized by iDLl but for molecular clouds this ~10 /im feature is too 
broad ( Ivan Breeman et al.ll201ll ). The competing hypotheses explaining the molecular cloud 
spectra could be tested by comparing with the glass spectra presented herein. 
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The importance of dust to astrophysical processes cannot be over stated. For example 
observations of high redshift {Z > 7) galaxies and quasars demonstrate tha t there was copi- 



ous dust produced by th e time the Universe was ~700 million years old (e.g.. lSugerman et al. 



20061 ; iDwek et al.l 120071 ). Furthermore, in Active Galactic Nuclei (AGN), the observed sili- 



cate absorption feature is shifted to peak at a longer wavelength, and is broader than that 
observed in o ur own Galaxy. T his spectral shift has been attribute d to calcium-aluminum 
rich silicates ( Uaffe et al.l 120041 ) or porous particles (ILi et al.l 120081 ). Moreover, data from 
Spitzer (SAGE-IRS Legacy program) show a number of both AGB and YSO sources with 
remarkably blue silicate features. Given that Mg, Fe and Si are formed in different nucle- 
osynthetic processes, the abundance ratios of these elements do not necessarily scale with 
metallicity. Disentangling dust formation mechanisms through observations of dust requires 
the optical properties of a range of silicate samples of varying Mg- Fe- and Si- contents and 
other components as provided here. The work presented here is a subset of a larger study 
and further sample compositions and structures will be presented in the near future. 



5. Conclusions 

We have presented new laboratory spectra of astrophysically relevant silicate glasses and 
compared them to existing data in the literature. We have shown that (1) "disordered" is 
not synonymous with glassy. In addition to structural disorder, porosity also affects spectral 
features. (2) Sample preparation and characterization are important. (3) We confirm the 
general trend of decreasing peak wavelength with increasing polymerization for the ~10 /im 
feature. However the scatter about this overall trend indicate that other compositional 
factors must be important. (4) Nothing quite matches the diffuse ISM in peak position and 
breadth. 

Spectral parameters of disordered silicates are sensitive to composition and sample syn- 
thesis techniques, which reflect degree of disorder, porosity, oxidation states and water con- 
tent. To understand dust formation we must disentangle these parameters through further 
systematic study of major compositional series using high resolution spectroscopy on thor- 
oughly characterized samples. Further studies of the parameters will follow. These new data 
can be used for interpretation of more esoteric environments e.g. high redshift galaxies and 
novae. 
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Fig. 1. — (a) Cartoon of possible grain structures, including single and polycrystalline grains, 
(b) Schematic of grain porosity vs structural disorder for different experimental sample ma- 
terials, (c) Simple phase diagram illustrating the location of the glass transition range, 
below the melting point of crystalline materials. Which side of the phase diagram is relevant 
depends on the cooling rate of the system. 
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Fig. 2. — Schematic structure of silicate glasses showing bridging oxygens joining tetrahedral 
units, and non-bridging oxygens coordinated by network-modifying cations. Trivalent tetra- 
hedral cations (Al 3+ , Fe 3+ ) must be charge-balanced by metal cations occupying interstices 
in the structure (e.g. Na + , Ca 2+ ~) The real three-dimensional structure of glass is a modified 
random network (e.g. iHendersonl 120051 ) . 



-30 - 



[1 760-1 300K] 
Al 2 3 , Alumina 



[1 625-1 400K] 
Ca 2 AI 2 Si0 7 , Gehlenite 



[1 450-1 100K] 
CaMgSi 2 6 , Diopside 



[1 360-1 000K] 
CaAI 2 Si 2 8 , Anorthite 



[1 440-1 050K] 
Mg 2 Si0 4 , Forsterite 



[1 350-1 040K] 
MgSi0 3 , Enstatite 



[1100-950K] 
[Mg,Fe] 2 Si0 4 , Olivine 



Fig. 3. — Predicted condensation sequence for O-rich environments jGrossmanlll972l l 
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Fig. 4. — New laboratory absorbance spectra of glasses. In all cases, x-axis is wavelength in 
jum; y-axis is absorbance. The dotted line indicates the barycentric position for each "10 /im" 
feature. The wavelength of these positions, along with the peak position and FWHMa are 
listed in Table HI The precise compositions for each sample are listed in Table EJ 
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Fig. 5. — New laboratory absorbance spectra where data extends to far-IR. In all cases, 
x-axis is wavelength in /im; y-axis is absorbance. The dotted line indicates the barycentric 
position for the "10 /im" and "18 /im" features. The wavelength of these positions, along 
with the peak position and FWHMa are listed in Table HI The precise compositions for each 
sample are listed in Table EJ Interference wide fringes are indicated by arrows for gehlenite 
remelt and enstatite, while narrow fringes are marked on the spectrum of Cosmic Silicate. 
The attributions for the observed features at ~ 10, 14, and 18/xm are indicated in the lower 
right panel. 
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Fig. 6. — Calculated Q values (absorption efficiency factors) for our samples In all cases, 
x-axis is wavelength in /xm; y-axis is Q a bs- The spectra are plotted such that the barycentric 
position of the feature is reddest at the top and bluest at the bottom and are offset for 
clarity. The values for the offsets in y, along with the barycentric positions, peak position 
and FWHMa are listed in Table H 
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Fig. 8. — Laboratory spectra demonstrating changes in peak p osition and s hape with 
composition, a) top p anel: the effect of SiC>2 content, data from Ijager et al.l (119941 ) and 



Dorschner et al.l (119951 ); b) middle panel: the effect of SiC>2 content, our data on Mg- 
endmember glasses; c) bottom panel: the effect of Fe content on "pyroxene" glasses, where 
En# = 100 x Mg/(Mg+Fe) in the glass. The iron-bearing glasses are substantially oxidized. 
Data from iDorschner et al.l (119951 ). 



10 15 20 



wavelength (/im) 

Fig. 9. — Comparison of "dirty" silicates. Solid black line is our "Cosmic" silicate 
[Fe/(Mg+Fe)=0]; dotted black line is our "Basalt" [Fe/(Mg+Fe)=0.39]; dotted-dashed light 
grey line is Ensta t ite Re melt [Fe/(Mg+Fe)=0.07] dashed dark grey line is the "dirty" silicate 
from iJager et al.l ( 119941 ) 
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Fig. 10. — Comparison of synthetic spectra from IDLI and lOHMl with new laboratory spectra 
on akermanite. x-axis is wavelength in /im; y-axis is the imaginary part of the complex index 
of refraction (k). Left panels: synthetic akermanite; right panels: akermanite remelt; top row: 
comparison to cool oxygen-rich silicate from OHM; middle row: comparison to warm oxygen- 
deficient silicate from OHM; bottom row: comparison to DL astronomical silicate. Solid lines 
are the synthetic spectra; dotted lines are the new laboratory spectra. 
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Fig. 11. — Comparison of synthetic spectra from iDLl and lOHMI with new laboratory spectra 
on gehlenite. x-axis is wavelength in //m; y-axis is the imaginary part of the complex index 
of refraction (k). Left panels: synthetic gehlenite; right panels: gehlenite remelt; top row: 
comparison to cool oxygen-rich silicate from OHM; middle row: comparison to warm oxygen- 
deficient silicate from OHM; bottom row: comparison to DL astronomical silicate. 
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Fig. 12. — Comparison of synthetic spectra from IdlI and OHM with new laboratory spectra 
on forsterite and enstatite. x-axis is wavelength in fim; y-axis is the imaginary part of 
the complex index of refraction (k). Left panels: forsterite; right panels: enstatite; top 
row: comparison to cool oxygen-rich silicate from OHM; middle row: comparison to warm 
oxygen-deficient silicate from OHM; bottom row: comparison to DL astronomical silicate. 
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Fig. 14. — Comparison of synthetic spectra from iDLl and lOHMl with new laboratory spectra 
on "Cosmic Silicate" and Diopside. x-axis is wavelength in /im; y-axis is the imaginary part 
of the complex index of refraction (A;). Left panels: Cosmic silicate; right panels: diopside; top 
row: comparison to cool oxygen-rich silicate from OHM; middle row: comparison to warm 
oxygen-deficient silicate from OHM; bottom row: comparison to DL astronomical silicate. 
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Table 1: Laboratory spectral data from astronomical literature 



Citation 



Sample 

Preparation 

Technique 



Sample Barycenter FWHM 
Composition (/im) (A* m ) 



Barycenter FWHM 
(H (^m) 



Davl (119791) 



chemical vapor 
condensation 



Scott fc Dulevl laser ablation of 
( 19961 ) crystalline samples 



Kratschmer fc Hcufferraapiation of 
(|l979l ) crystalline samples 



Jena 



Melting and 
Quenching 



Dorschner et aLI Melting and 
|l995h Quenching 



FolOO 
EnlOO 
FolOO 
EnlOO 
FolOO 

Fo50 

En50 
EnlOO 



10.2 
10.0 
10.4 
10.0 
10.5 

10.3 

9.8 
9.8 



Jager et al 
( 20031) 



Jager 



(1994) 



et al 



(119981) 



DL 



OHM 



OHM cold 



1.9 
1.8 
1.8 
1.8 
2.2 

2.3 

2.6 
1.8 



21.8 
20.9 
22.3 
22.5 
18.0 

20.5 

20.0 
20.4 



11.2 

6.9 

17.7 

17.3 

2.5 

11.5 

7.3 
9.8 





En80 


9.8 


2.3 


20.2 


9.1 




En60 


9.9 


2.5 


20.3 


9.3 




En40 


9.7 


2.7 


20.4 


8.2 


Sol-Gel 


FolOO 


10.0 


2.2 


17.7 


4.0 




EnlOO 


9.5 


2.0 


20.3 


9.1 


Melting and 


dirtysil 1 


9.6 


1.4 


21.3 


10.2 


Quenching 












Melting and 


Gehlcnitc 


10.2 


1.2 


21.4 


7.5 


Quenching 














mclilitc 


10.5 


1.8 


19.9 


5.6 


synthetic 


generic 


10.5 


3.3 


20.8 


8.6 


synthetic 


generic 


10.7 


2.5 


19.7 


5.7 


synthetic 


generic 


10.8 


2.6 


19.7 


5.3 



1 Mg .5oFeo.43Cao.o3Alo.o4Si0 3 
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Table 2: Glass Sample Compositions and physical parameters 



Name 


Formula 


P (kg m~ 


3 ) T g (K) A1/(A1+ 


Si) NBO/T Q 


Synthetic Gehlenite 


(Nao.oiCai.99)Al(Ali.o3Sio.98)0 7 


2879 


1111 


0.68 


0.65 


Gehlcnite Remelt 


Nao.02Ca1.97Mgo.37Feo.06Al1.14Si1.44O7 


2927 


1032 


0.44 


1.43 


Synthetic Akermanite Ca2.04Mgo.9sAlo.02Si1.97O7 


2951 


1012 


0.01 


3.04 


Forsterite 


Mg2.035Sio.98304 


2920 


1040 fc 


0.00 


4.00 


Diopside 


Ca1.02Mgo.92Si2.03O6 


2853 


999 


0.00 


1.91 


Akermanite Remelt 


Na0.30Ca1.70Mg0.47Fe0.29Al0.54Si1.64O7 


2995 


932 


0.25 


2.15 


Basalt 


K0.0iNa0.06Mg0.45Ca0.54Fe0.32Ti0.09Al0.70Sil.72O6 






0.03 


1.48 


Synthetic Enstatite 


(Mgi. 98 Alo.o2)(Alo.oiSii.99)0 6 


2735 


1037 


0.02 


1.94 


Enstatite remelt 


Ca0.02Mg1.8iFe0.14Al0.03Si2.00O6 


2769 


1022 


0.01 


1.93 


"Cosmic Silicate" 


(Nao.iiCao.i2Mgi. 86 )(Alo.i8Sii.85)0 6 


2772 


1001 


0.09 


1.92 


Hcrasil 


Si0 2 [> 0.1wt%H 2 O] 


2200 c 


1420 


0.00 


0.00 


Obsidian 


Ko.28Nao.29Cao.02Feo.04Alo.64Si3.35O8 






0.16 


0.02 



a NBO/T = (K + Na + 2Ca + 2Mg + 2Fc 2 + - Al - Fc 3 +)/(Si + Ti + Al + Fc 3 +) 
b from Tangcman ct al. 



c from manufacturer 
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Table 3: Compositions of glasses, determined by electron microprobe, in weight percent 
oxides 



name 


Si0 2 


TiOz 


A1 2 3 


weight percent oxide 

FeO MnO MgO CaO Na 2 K 2 P 2 5 


total 


Water content" 
(ppm) 


Synthetic gchlcnitc 


21.53 


na 


37.88 


na 


na 


na 


40.85 


0.12 


na 


na 


100.37 


76 


Gchlenite remelt 


31.46 


0.01 


21.11 


1.58 


0.01 


5.47 


40.26 


0.27 


0.02 


0.04 


100.22 


85 


Synthetic akermanite 


43.12 


na 


0.44 


na 


na 


14.33 


41.62 


na 


na 


na 


99.51 


80 


Forsterite 


41.17 


na 


na 


na 


na 


57.19 


na 


na 


na 


na 


98.36 


bd 


Diopside 


54.82 


na 


na 


na 


na 


16.55 


25.58 


na 


na 


na 


96.95 


57 


Akermanite remelt 


34.34 


0.89 


9.65 


6.85 


0.52 


6.53 


33.32 


3.08 


0.49 


2.03 


97.72 


125 


Basalt 


45.92 


3.04 


15.99 


10.16 


0.13 


8.04 


13.45 


0.80 


0.18 


na 


97.71 


bd 


Synthetic enstatite 


59.69 


na 


0.68 


na 


na 


39.88 


na 


na 


na 


na 


100.25 


80 


Enstatite remelt 


58.30 


0.02 


0.66 


4.74 


0.12 


35.43 


0.54 


0.02 


0.01 


na 


99.85 


63 


Cosmic silicate 


54.26 


na 


4.34 


na 


na 


36.58 


3.27 


1.61 


na 


na 


100.05 


81 


Obsidian 


74.82 


0.03 


12.09 


0.95 


0.03 


0.01 


0.39 


3.39 


4.91 


na 


96.62 


-3000 


Hcrasil fc 


100.00 


na 


na 


na 


na 


na 


na 


na 


na 


na 


100.00 


940 



na = not analyzed (mostly for synthetic melts where the ingredients were well known) 
bd = below detection limit 

a determined from near-IR spectra using the method in Hofmcister et al. (2009) 
6 pure Si0 2 with 15-26 ppm metals (Hofmeister and Whittington, in review) 
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Table 4: Spectral parameters of sample glasses. In all cases the parameters are measured 
from the absorbance data and the Q^ s data as plotted in Fig. 0] and Fig. El respectively 



SAMPLE 




~10 /im 








~18 /im 






10/18/zm 




NAME 


Peak 


Barycenter 


FHWM 


Peak 


Barycenter 


FHWM 


strenj 


;th ratio* 


Offset 1 ' 




( H 


(/im) 


a 


Q 


(H 


(H 


a 


Q 


a 


Q 




Gehlcnite 


10.3 


10.8 


2.71 


2.86 














+3.25 


Gehlenite RM 


10.8 


10.6 


2.69 


2.88 


_} 


20.0 


4.08 


4.45 


3.74 


4.25 


+2.75 


Akermanite 


10.3 


10.6 


2.31 


2.46 














+2.50 


Forsterite 


10.2 


10.4 


2.25 


2.43 














+2.25 


Diopside 


10.1 


10.2 


2.33 


2.46 


19.2 


19.4 


4.34 


4.12 


1.65 


1.26 


+ 1.50 


Akermanite RM 


10.3 


10.1 


2.63 


2.70 


21.0 


20.9 


4.45 


4.70 


1.78 


1.57 


+ 1.25 


Basalt 


10.0 


9.9 


2.27 


2.41 














+ 1.00 


Enstatite 


9.9 


9.9 


2.38 


2.58 


17.6 


18.5 


3.99 


4.70 


3.03 


2.25 


+0.75 


Entatite remelt 


10.0 


9.8 


2.33 


2.48 














+0.50 


Cosmic silicate 


9.8 


9.7 


2.76 


3.15 


18.3 


18.7 


4.64 


4.54 


4.87 


4.79 


+0.25 


Obsidian 


9.0 


9.1 


1.91 


2.14 














0.00 


Herasil 


9.0 


8.9 


1.24 


1.34 














-0.5 



* strength is denned as the equivalent width (EW) of the feature, rather than the peak-to-continuum ratio. 
Therefore the strength ratio is EWio^m/ EWis Mm - 

Peak position for the Gehlenite remelt sample cannot be accurately measured because of interference fringes, 
t the offset is added to the normalized Q- values and is applied in Figure [5] for clarity, 
t strong fringing prevents measurement of peak position. 
RM designates remelted samples (see § 12. 2} 



